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Resonant second harmonic generation between 1550 nm and 775 nrn with normalized outside efficiency > 
3.8 x 10~ 4 mW -1 is demonstrated in a gallium phosphide microdisk supporting high-Q modes at visible 
(Q ~ 10 4 ) and infrared (Q ~ 10 5 ) wavelengths. The double resonance condition is satisfied for a specific 
pump power through intracavity photothermal temperature tuning using ~ 360 ffW of 1550 nm light input 
to a fiber taper and coupled to a microdisk resonance. Power dependent efficiency consistent with a simple 
model for thermal tuning of the double resonance condition is observed. 


Since the first observation of second harmonic gen¬ 
eration (SHG) in 1961 1 , it has become a ubiquitous 
demonstration of nonlinear optics. In recent years there 
has been mounting interest in SHG within micron-scale 
optical structures such as waveguides and cavities 2-12 . 
These experiments seek to take advantage of devices 
whose combination of large optical quality factor, Q , 
and small mode volume, F, provide enhancements to 
electromagnetic per-photon field intensities at both fun¬ 
damental and second harmonic wavelengths. Further¬ 
more, the compact nature of these optical devices lends 
itself to convenient integration into complex on-chip pho¬ 
tonic circuits. To date SHG has been demonstrated 
in a number of microresonator geometries including 
microdisks 7,11-13 , microrings 14 , microspheres 15 , photonic 
crystal nanocavities 5,6,8,9 and waveguides 16,17 . An im¬ 
pressive 9x 10 -2 mW _1 conversion efficiency was achieved 
by Furst et. al. in macroscopic (mm) sized whispering 
gallery mode resonators 18 . Conversion between 1985 nm 
and 993 nm light with normalized outside and circulating 
efficiencies of 10~ 5 mW -1 and 10 _3 mW -1 , respectively, 
have been demonstrated in /im-scale gallium arsenide mi¬ 
crodisks with low optical absorption at IR wavelengths 12 . 
Microcavities have also been utilized in nonlinear opti¬ 
cal processes including sum-frequency generation 19 , and 
third harmonic generation 20 . 

In this work we study microdisk cavities such as 
the device shown in Fig. 1(a), fabricated from gallium 
phosphide (GaP). The optical transparency window of 
GaP spans wavelengths from 550 nm to IR, making it 
a promising material for nonlinear wavelength conver¬ 
sion between C-band (1550 nm) and visible wavelengths 8 . 
Second harmonic conversion efficiency can be enhanced 
by fabricating microcavities supporting high-Q optical 
resonances at A c and A c that are resonant with the pump 
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FIG. 1. (a) SEM of a GaP microdisk before being undercut, 
(b) generated by the m = 27 fundamental TE mode at 
A c ~ 1550nm of a microdisk with dimension as in (a). 


(A 0 ) and second harmonic wavelengths (A c = A 0 /2), re¬ 
spectively. The high sensitivity of small-F devices to fab¬ 
rication imperfections and variations in material optical 
properties, combined with their large free-spectral range, 
makes it challenging in practice to realize high-Q/F 
doubly-resonant microcavities. Approaches for address¬ 
ing this problem include near-field perturbative tuning 21 , 
and bulk temperature tuning 13 . Here we demonstrate a 
method for tuning optical resonance wavelengths in-situ 
using dispersive intracavity thermo-optic effects, and uti¬ 
lize this tuning to enable highly-efficient conversion be¬ 
tween A 0 = 1545 nm and A 0 = 772 nm, with normalized 
outside efficiency in excess of 3.8 x 10~ 4 mW -1 . 

Efficient SHG in microcavities requires, in addition 
to resonances at A 0 and A 0 , a phasematching mecha¬ 
nism to overcome dispersion intrinsic to the microcav¬ 
ity. A common approach to achieve this, known as 
quasi-phasematching, relies upon periodic domains with 
alternating nonlinear susceptibility 22 . The zincblende 
structure of GaP possess 4 symmetry 22 , which can be 
used to realize quasi-phasematcliing of microdisk whis¬ 
pering gallery modes without explicit creation of periodic 
domains 10,12,23 . The microdisks studied here support ra¬ 
dially polarized (TE) modes near A 0 , and both TE and 
vertically (z) polarized (TM) modes near A OJ that are 
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coupled by the second order nonlinear susceptibility of 
GaP. The amplitude of the radial and z field components 
these modes varies with e lmS , where m is the azimuthal 
mode number and 9 is the cylindrical angular coordi¬ 
nate. The corresponding cartesian in-plane components 
E x and E y of the TE modes experiences a sign change 
with period |A0| = n. As a result, the second order non¬ 
linear polarization along z, L = 2eodi^E x E y where di 4 
is the relevant nonlinear susceptibility tensor element of 
GaP, experiences a sign change with period |A0| = tt/ 2. 
This is visualized in Fig. 1(b), and can be interpreted 
as momentum imparted by the periodic effective inver¬ 
sion of crystal orientation, creating azimuthal momentum 
components at 2 to ±2. As a result, the time-averaged 
coupling between P 1 ) 11 - and the TM microdisk mode with 
azimuthal index to is maximized when to = 2m±2. This 
±2 offset can compensate for microdisk dispersion 10 , so 
that satisfying \ c \ m = 2 A c | 2 m ±2 becomes possible. 

The impact of these effects on SHG is captured by 10,24 : 

(1) 

2 

p\ 

which describes the second harmonic power P output into 
a waveguide coupled to the microdisk for IR pump power 
P input to the waveguide. Q e j and Q Cj t are the external 
waveguide coupling (e) and the total (t) quality factors 
of the pump and IR microdisk modes, respectively. The 
quasi-phasematching is captured by the second harmonic 
coefficient, K, given in the Supplementary Material 24 . 
The remaining terms describe the cavity enhancement, 
and are maximized at “double resonance” where both the 
pump and second harmonic fields are resonant with a cav¬ 
ity mode, i.e. A =^A = 0. Here, A(A 0 , Pd) = A 0 - A c (P d ) 
and A(A 0 ,P d ) = A 0 - A c (P d ) = A 0 /2-A c (P d ) are the de¬ 
tunings between the pump and second harmonic fields, 
respectively, and the microdisk modes. Due to pho- 
tothermal effects discussed below, A c and A c depend on 
the pump power dropped into the cavity, P d . 

In microdisks, double resonance is generally not sat¬ 
isfied intrinsically. However, adjusting the microdisk 
temperature, T, can compensate for an initially non¬ 
zero relative detuning |A C — 2A C |. Changes in T modify 
A c and A c via thermal expansion and the thermo-optic 
effect 25 , and tuning rates d\ c /dT and d\ c /dT can differ 
due to dispersion of the modal confinement, refractive in¬ 
dex (neap)) and thermo-optic coefficient. In the experi¬ 
ment described below, the wavelength dependence of the 
normalized thermo-optic coefficient (l/nQ a p)dnQ a p/dT 
is the dominant mechanism for tuning rate dispersion. 
For GaP we calculate this coefficient at A c and A c to be 
3.4 x 10" 5 [K” 1 ] and 2.9 x 10" 5 [K” 1 ] respectively 26 , 
which leads to a differential tuning coefficient r] = 1 — 
(2d\ c /dT)/{d\ c /dT) = 0.176. 
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FIG. 2. (a) Cartoon of influence of P on detunings A(A D ; P ) 
(blue) and 2A(A 0 ;P) (red) for rj > 0. The dashed lines show 
the expected behavior in absence of thermo-optic effects, (b) 
IR transmission spectrum of a fiber taper coupled to a 6.52 pm 
diameter microdisk, (c) Red line: Visible transmission spec¬ 
trum of the fiber taper. Green line: SHG signal generated 
by the IR input light from (b) and collected by the fiber 
taper. Note that for the P used in (b) and (c), the short¬ 
est wavelength IR mode nearly satisfies the double resonance 
condition. 


In our experiment we do not have independent control 
of T. However, because of linear optical absorption and 
accompanying heating, T is proportional to P d ( A). For a 
given 77 and initial |A C —2A C |, it is possible to find a critical 
power P = P cr it where double resonance is achieved. The 
principle of this scheme is illustrated in Fig. 2(a), where 
we sketch A and 2A versus A 0 for three values of P, 
assuming ij > 0. As A 0 is tuned towards A c from blue 
to red, the increase in P d { A) heats the cavity, causing 
both A c and A c to shift to longer wavelengths. For P < 
Pcrit (left panel in Fig. 2 (a)) the deviation is relatively 
small and double resonance is never realized. When P = 
-Pcrit (center panel in Fig. 2 (a)), the double resonance 
condition is satisfied. For P > P cr j t (right panel in Fig. 
2 (a)), A c is shifted past 2A C before A 0 reaches A c , and 
double resonance is never achieved. A sharp change in 
A and A occurs in all three scenarios when A 0 “catches 
up” to the thermally shifted A c 25 . 

We now show that this tuning scheme can be used 
to demonstrate SHG at double resonance in GaP mi¬ 
crodisks fabricated following Mitchell et al. 2 '. To iden¬ 
tify promising microdisks for efficient SHG, devices were 
characterized using fiber taper mode spectroscopy at IR 
wavelengths, while the spectrum of SHG produced by 
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FIG. 3. (a) Normalized transmission of the pump through 
the fiber taper when A 0 is scanned from 1544 nm to 1548 
nm (in the direction of increasing A 0 ) for various fixed values 
of P. A thermally induced cavity resonance shift is evident 
with increased P. The red dashed line traces over points 
of minimum transmission where A 0 = A c . (b) SHG signal 
corresponding to the operating conditions in (a). The dashed 
line traces A c from (a), and serves as a guide to the eye. 


the microdisk and collected by the fiber taper was si¬ 
multaneously measured. Figure 2(b) shows the trans¬ 
mission through the fiber taper when it is positioned in 
the near field of a 6.52 pm diameter microdisk, measured 
with P ~ 0.5 mW from a tunable IR wavelength laser 
(New Focus TLB-6700). Evanescent coupling from the 
fiber taper to high-Q microdisk doublet modes 28 results 
in sharp transmission dips and corresponding large Pd- 
Through comparison with finite difference time domain 
(MEEP) simulations of the microdisk mode spectrum 
measured from 1470—1570 nm, we identify doublets near 
1545 nm and 1554 nm as corresponding to TE-polarized 
modes with m, radial (p), and axial ( q ) numbers {27,1,1} 
and {19, 3,1}, respectively. From fits to the doublet line- 
shape, we measured unloaded Q of 1.1 x 10 5 and 2.8 x 10 4 
for the p = 1 and p = 3 doublets, respectively. Figure 
2(c) shows that SHG is observed when pumping both 
modes, and is strongest for the p = 1 mode. 

To study the initial alignment between modes near 
A c and the SHG signal, we measured the fiber taper 
transmission using a visible supercontinuum source and 
spectrometer detection. Comparing the transmission and 
SHG spectra in Fig. 2(c) reveals that the strongest SHG 
occurs when exciting a doublet mode with SHG emission 
at A 0 = A 0 /2 that is close to a resonance at A c ~ 772.2 nm 
({m,p, q} = {56,3,1}, Q ~ 9600). In contrast, SHG 
from the other mode in the doublet is weak, as it is off 
resonance from A c . 


Further insight into the role played by double reso¬ 
nance was provided using the tuning scheme described 
above. We scanned A 0 from blue to red across A c for 
a range of fixed values of P and measured the pump 
transmission (Fig 3(a)), and the corresponding SHG (Fig. 
3(b)) as a function of P and A 0 . Each point in Figs. 3(a) 
and 3(b) is the integrated intensity of the narrowband 
signal detected by the IR and visible spectrometers at A 0 
and A 0 respectively, for a given {A C ,P}. The redshift of 
A c in the pump transmission spectrum in Fig. 3(a) with 
increasing P results from thermo-optic effects in the mi¬ 
crodisk. These effects are also responsible for the asym¬ 
metric “shark fin” lineshape evident in Fig. 2(b) 25 . Note 
that the jagged edges in Fig. 3(a) for high P are due to 
the discrete stepping of P. Figure 3(b) shows that SHG is 
most intense when exciting the shorter wavelength mode 
of the doublet, consistent with the observation in Fig. 
2(c). Additionally, SHG is maximized at an intermedi¬ 
ate value of P. As described below, this can be explained 
by the Pd dependence of A(A 0 = A c ) illustrated in Fig. 
2(a). 

The SHG behaviour can be investigated more quanti¬ 
tatively by analyzing the subset of data in Fig. 3 where 
the pump is on-resonance (A = 0). In this case the SHG 
signal is described by a simplified form of Eq. (I) 10,24 : 


P(A = 0) oc 


(Ac/2 Qt] 2 


A(A 0 ,P) 2 + (A c /2Q t )' 


x |A| 2 P 2 (2) 


This expression is a product of two functions: the first 
captures the impact of the P dependence of the detuning 
A between the SHG signal (A 0 = A c /2) and the nearest 
cavity mode (A c ), and the second describes the usual P 2 
dependence of SHG, and scales with the effective non¬ 
linear susceptibility of the microdisk modes of interest. 
Conveniently, A is a linear function of P and A 0 for A = 0 
(see Supplementary Information 24 ), allowing Eq. (2) to 
be easily fit to experimental data. 

In Fig. 4 we analyze the A = 0 data from Fig. 3 us¬ 
ing Eq. (2). We only consider the shorter A c doublet 
mode, and find A C (P; A = 0) from the data in Fig. 3(a). 
In Fig. 4(a), we plot SHG absolute efficiency vs. P for 
A = 0. Here absolute efficiency is defined as P/P, where 
powers are measured in the fiber taper immediately be¬ 
fore (P) and after (P) the microdisk, taking into account 
asymmetric fiber taper insertion loss. We emphasize that 
because of the P dependence of A c shown in Fig. 3(b), 
for every value of P, the selected A 0 corresponding to 
A = 0 varies. For low P < 0.18 mW, the absolute effi¬ 
ciency is observed to increase approximately linearly with 
P, as thermo-optic effects are small compared to the in¬ 
trinsic “cold cavity” A|p = o- In this region we estimate 
A to range from —34 pm to —28 pm. As P approaches 
P = Pcrit = 0.35 mW, thermo-optic effects become sig¬ 
nificant, A —>■ 0 satisfying double resonance, and the 
absolute efficiency increases superlinearly with P to a 
maximum value of 1.5 x 10 -4 . For P > P cr it, absolute 
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FIG. 4. (a) SHG absolute efficiency P/P versus P for A = 0. 
The red line is a fit to the data using Eq. (2). The inset 
highlights low P data and includes a linear fit. (b) Normalized 
SHG count rate (oc P/P 2 ) plotted versus pump wavelength 
A„ for A = 0, where the red line is a fit to the data adapted 
from Eq. (2). (c) SH power vs. P for A = 0. The green points 
show the unprocessed measured data. The blue points show 
the predicted SH power if A = 0 for the entire range of P. 
The red line is a weighted least squares fit to the data of a 
function proportional to P 2 , and indicates a normalized SHG 
efficiency P/P 2 = 3.8 ± 0.2 x 10 -4 mW -1 . 


efficiency is observed to decrease. 

The observed P dependence of P(A = 0) is consistent 
with A (P; A = 0) increasing monotonically with P from 
A < 0 to A > 0. We test this quantitatively by com¬ 
paring the predictions from Eq. (2) with the measured 
data in Fig. 4(a). Here the required fitting parameters 
are an overall scaling factor, and the differential thermo¬ 
optic tuning coefficient 77 that determines the linear de¬ 
pendence of A on P for A = 0. We find a best fit for 
?y = 0.181, which is within 3% of the theoretical value 
presented above. In applying this model, we included a 
A 0 dependence in the scaling factor that follows an Airy 
function and accounts for etaloning in the non-wedged 
neutral density filters used in this experiment, resulting 
in oscillations evident in the data. 

Further analysis of the SHG signal at A = 0 is given 
in Fig. 4(b), which shows the SHG signal for varying A 0 . 
Here we are plotting P normalized by P 2 in order to 
isolate the cavity contribution described in Eq. (2). The 


resulting fit of Eq. (2) to the data, shown in Fig. 4(b), 
has good agreement. For this fit, we use the value for 
?y found above, and treat Q t as a fitting parameter; the 
fit predicts Qt ~ 1.0 x 10 4 , in close agreement with the 
spectrometer measured value in Fig. 2(c). As in the P- 
dependent fit in Fig. 4(a), we included a A 0 dependent 
scaling factor to account for etaloning from the neutral 
density filters that creates the observed oscillations. 

Finally, in Fig. 4(c) we analyze the SHG signal power 
vs. P at A = 0. In this analysis, we plot P in absolute 
measured units, as well as scaled by the P dependent 
microdisk resonance response described by the first term 
in Eq. (2) and predicted from the fit to A(P; A = 0) ob¬ 
tained for Fig. 4(a). This latter scaling isolates the term 
in Eq. (2) that describes the microdisk nonlinear suscep¬ 
tibility, and is ideally proportional to \K\ 2 P 2 . As shown 
in Fig. 4(c), when scaled in this manner, the data is well 
described by a P 2 dependence. This plot represents the 
predicted SHG signal if the enhancement provided by the 
microdisk density of states is fixed to its maximum value 
at A = A = 0. Note that in our experiment, this ef¬ 
ficiency is only realized when P = P CI it, as evident by 
comparing the unsealed and predicted data in Fig. 4(c). 

The maximum normalized outside efficiency, defined as 
P/P 2 , was found from the point of highest efficiency in 
Fig. 4(a) to be 4.4 x 10~ 4 mW -1 for the case P = P C1 ; t . 
The weighted least squares fit to the scaled data shown 
in Fig. 4(c) can be used as a consistency check, and gives 
the predicted normalized outside efficiency if the dou¬ 
ble resonance condition A = A = 0 is satisfied over 
the entire range of P. We find a normalized efficiency 
of (3.8 ± 0.2) x 10 -4 mW -1 , in good agreement with the 
measured maximum normalized outside efficiency at P crit 
in Fig. 4(b) where the double resonance condition is sat¬ 
isfied. In comparison, the SHG signal at low P shown 
in the inset Fig. 4(b), where A is approximately con¬ 
stant (i.e. independent of P), has a normalized outside 
efficiency of 1.78 x 10 -4 mW -1 . This illustrates that the 
SHG normalized outside efficiency P/P 2 , in addition to 
the total efficiency P/P shown in Fig. 4(b), is enhanced 
through satisfaction of the double resonance condition at 

P = Petit- 

In conclusion, we have demonstrated resonant sec¬ 
ond harmonic generation from 1550 nm to 775 nm in 
a GaP microdisk with maximum normalized outside ef¬ 
ficiency of 3.8 x 10 -4 mW -1 for P = P cr jt = 0.35 mW, 
which is larger than previously reported values in simi¬ 
larly sized structures used for 2000 nm to 1000 nm wave¬ 
length conversion 12 . We have shown that this efficiency 
is achieved via double resonance between high-Q modes 
of the GaP microdisk at both IR (Q ~ 1.1 x 10 5 ) and vis¬ 
ible (Q ~ 1.0 x 10 4 ) wavelengths. Further improvements 
to the efficiency could be made by spectrally aligning 
higher-Q modes, and optimizing fiber taper coupling at 
visible wavelengths. Control of the differential tuning in¬ 
dependent of pump power could be realized by thermally 
heating using light from an additional optical mode not 
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involved in the nonlinear conversion process. As GaP is 
a piezoelectric material, it may be possible to tune the 
modes through electronic means 29 . The geometry pre¬ 
sented in this letter can also be adapted for use in other 
nonlinear optics scenarios such as four wave mixing and 
downconversion 30-33 . With the present system we have 
already observed third harmonic generation, although ef¬ 
ficient collection of the resultant 515 nm light via fiber 
taper requires additional optimization. 
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I. THERMAL TUNING 

The thermal dependance of the laser cavity detuning 
can be approximated by the expression 25 : 

A(ST) « A 0 — A c [1 + aST ], (SI) 

with a = e + 4 ^. Here ST is the deviance from the 
equilibrium temperature, e is the thermal expansion co¬ 
efficient, and n is the refractive index. 

Furthermore, we can write an expression for the change 
in cavity temperature due to optical absorption when the 
cavity is pumped at resonance and is in thermal equilib- 


AT = 


2 Q 2 Pf p 

QnQabsK th tk h 


where Q a bs is the quality factor associated with opti¬ 
cal absorption in the cavity, and K t h is the inverse of 
the thermal time constant of the cavity. This expres¬ 
sion is derived by considering the rate at which power 
is absorbed into the cavity when optically pumped on 
resonance, and equating it to KthST, the rate at which 
power is dissipated. For simplicity we have grouped to¬ 
gether terms into the constant Cth 

Making use of equations (S2) and (SI) we can write 
expressions for A for the case when A = 0: 


a \\ 0 + ciq 
biPf + b 0 , 


where: 


clq — |^1 A 0 , 

1 a\ 0 

Ql= 

, Mo ~~ 

— Cth 2 CtA 0 , 


bi = Y~ x °- 


II. SHG IN DOUBLETS 

The optical resonance used in the current experiment 
was a doublet mode. Such modes occur in whispering- 
gallery type resonators when surface roughness causes 
the coupling rate between clockwise and counter clock¬ 
wise propagating mode to exceed the combined loss rate 
to all other channels 28 . The transmission spectra of the 
doublet contains two distinct resonances, each of which 
is a standing wave created from linear combinations of 
clockwise and counter clockwise propagating modes. In 
the case of the doublet considered in the current experi¬ 
ment, these standing waves are spectrally well separated, 
which permits the individual excitation of each of these 


modes. We write the amplitude of the standing wave at 
the fundamental frequency cug, as 28 : 


—k/ y/2S f 

i{A u ) - 7t/2 ’ 


(S 8 ) 


where A is the amplitude of the standing wave mode, Sf 
is the amplitude of the forward propagating mode in the 
fiber taper, k is the waveguide-cavity coupling, A w is the 
laser-cavity detuning in terms of frequency, and 7 1 is the 
energy decay rate of the cavity into all channels. Here 
|H | 2 is normalized to energy, and |<S /| 2 is normalized to 
power. 

In contrast to the travelling wave fundamental modes 
considered elsewhere 10 , a standing wave fundamental 
mode will excite both clockwise and counter clockwise 
modes at the second harmonic frequency tu 0 = 2 lu 0 . The 
coupled mode equations describing this process may be 
written as: 


dn cw ~ 

— = (iu> c — j t /2)a cw + C, (S9) 

at 

da ^ = {iu c -l t /2)a ccw + 1 (S10) 

at 


where £ is the second harmonic source term, 7 : t is the 
energy decay rate of the loaded cavity around the cavity 
frequency u> c , and a cw ’ ccw = A cw ' ccw e lu ° t with w 0 repre¬ 
senting the laser frequency and A cw,ccw as the amplitude 
of the clockwise or counterclockwise propagating modes. 
Note that in our experiment, the second harmonic mode 
was found to be a singlet mode, so surface roughness did 
not appeciably couple the clockwise and counter clock¬ 
wise propagating modes at w c . 

The explicit form of the nonlinear source term may be 
written as: 

(Sll) 


£ = 2-7T-— 

s 27 tR 


dA cw,ccw 

dd 


where v g is the group velocity of the SH mode, and R is 
the radius of the disk. The first term of this expression 
gives the frequency at which the SH mode circulates the 
disk in units of rad/s, and the second term describes the 
average gain per cycle. 

The nonlinear source term can be calculated from the 
spatial profile of the optical modes as: 


= A 2 5 u>fsr r A t( Am+2)8 + K _ e ~z(A m +2)0] dQ 
4tt 7 0 L J 

(S!2) 

where £ is the source term, Slursr is the f ree spectral 
range of the microdisk at ui 0 , and K + , K_ are the sec¬ 
ond harmonic coefficients, and A?n is the mismatch in 
the azimuthal numbers of the microdisk at w 0 and wo, 
defined as Am = m — 2m. The explicit form of the sec¬ 
ond harmonic coefficients, as derived in 10 , may be written 



2 



FIG. SI. (a) Second harmonic signal for the special case A = 0, as a function of Pf and A. The coloured curves represent 
possible trajectories of an experiment, (b) Power dependence of second harmonic generation for the case of constant detuning 
(magenta curve) and for variable detuning (green curve), (c) Detuning dependence of second harmonic generation of constant 
power (yellow curve), and for varied power (green curve). 


as: 


K+ = 


du r°° p ~ (ml 

2e 0 u} 0 n 4 J 0 J_ h/2 ^\ r 


ap) d - ip • 

(S13) 


du r r h/2 t (mi 
2e 0 uj 0 n 4 J 0 J- h / 2 ^ V r 
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dzdp , 
(S14) 


where n is the refractive index of GaP at ujq, and ip, ip give 
the radial and vertical dependence of the z-component 
of the inodes at wo and to o- Here we have assumed 
the normalization conditions f |?/>| 2 gL4 = Sujf sr /2TT and 
f |'i/’| 2 dH = Su}f sr /2j: where the integrals are taken over 
the infinite half-plane 10 . 

Finally we use S9, S10 to solve for the steady state 
amplitudes of the second harmonic modes propagating 
in the clockwise and counter clockwise directions: 


where is the laser-cavity detuning w 0 /2 — uj c . 

Making use of the amplitude normalizations, and 
waveguide-cavity coupling parameters, we can relate 
Pf.b-, the power in the taper at A 0 in either the forward 
or backward direction and Pf, the power in the taper at 


A 0 to the amplitudes: 


Pf = \Sf\ 2 , 

(S17) 

p ,2 W 2 \Af, b \ 2 

Pf,b - \£>f : b\ = 9 

(S18) 


Combining equations (S8), (S12), (S15), (S17), (S18) 
and writing the decay rates in terms of their respective 
quality factors we arrive at a final expression in terms of 
P/p and Pf. 


Pf,b = \K\ 2 


X C /2Q K 

A 2 + (A c /2 Q t y 


1 A e /2 Q k Y 2 

2 A 2 + (A c /2Qt) 2 J f 

(S19) 


j^ cw _2_ 

*A W - 7t/2 ’ 

j^CCW _C_ 

*A W - 7t/2 ’ 


(515) 

(516) 


with: 

\ K \ 2 = ^SX 2 FSR [K + e l ^ Am+2) smc[n(Am + 2)] 

+ AAe* 7r(Am - 2) sinc[7r(Am - 2)]] 2 . (S20) 
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In order to better visualize the physics behind the data 
analysis of the main text, we have plotted Eq. (S19) in 
Fig. SI for the case A = 0. In this instance, the expres¬ 
sion for Pf is seperable into a function of Pf and a func¬ 
tion of A. We display this on Fig. SI along with three 
curves representing various paths in Pf — A space which 
an experiment might take. The magenta curve represents 
the case of fixed detuning A = 0, which demonstrates the 
typical quadratic dependence on Pf expected from a sec¬ 
ond harmonic experiment. The yellow curve is the case 
of fixed power, with a variable detuning, which displays 
the Lorentzian line shape contributed by the cavity res¬ 
onance lineshape. Lastly, the green curve represents the 
case where both detuning and pump power are varied. 
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